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SUMMARY 

Realizing  the  effect  of  network  improvement  on  the 
operation  of  a  highway  transportation  system,  this  study 
develops  a  model  which  enables  us  to  obtain  an  optimum 
investment  policy  for  improving  a  transportation  network. 

Consider  an  area  with  the  future  traffic  demand  assumed 
to  be  known.   The  problem  then  is  to  build  (or  to  improve) 
a  transportation  network  to  accommodate  the  demand  with 
mimimum  overall  cost,  assuming  that  the  desired  network 
configuration  has  already  been  determined.   Travel  time 
cost  is  assumed  to  be  the  only  significant  factor  of 
operating  cost.   Therefore,  the  overall  cost,  or  trans- 
portation cost,  is  assumed  to  be  equal  to  the  sum  of  in- 
vestment cost  and  travel  time  cost. 

A  non-linear  total  travel  time  equation  was  developed 
to  express  the  travel  time  as  a  function  of  the  investment 
and  of  traffic  volume.   Three  formulations  associated  with 
three  different  investment  circumstances  are  presented. 
With  the  zonal  interchanges  considered  as  given  information, 
the  discrete  maximum  principle  (1)  was  utilized  to  assign 
the  trips  and  the  investment  to  each  link  of  the  network 
based  on  the  criterion  of  minimizing  the  overall  cost.   The 
maximum  principle  enables  the  location  of  the  optimum  with- 
out undue  difficulties. 


INTRODUCTION 

The  economic  analysis  of  a  transportation  network 
provides  valuable  guidance  in  developing  a  comprehensive, 
long-range  transportation  plan  which,  as  concluded  by 
Zettel  and  Carll,  is  the  basic  objective  of  a  transpor- 
tation study  (2).   Being  part  of  the  public  services  and 
competing  for  the  use  of  limited  resources,  the  transpor- 
tation system  should  be  built  and  operated  economically 
while  at  the  same  time  it  should  meet  the  standards  and 
goals  of  the  community  in  order  to  promote  growth  and 
meet  the  needs  of  the  economic  activities.   Specifically, 
the  objectives  of  a  transportation  system  have  been  sum- 
marized as:  (3) 

1.  Provide  a  means  for  moving  people  and  goods 
safely,  freely  and  economically. 

2.  Provide  a  choice  of  mode. 

3.  Make  the  city  a  more  attractive  place  to  live. 

4.  Provide  the  means  for  fulfilling  the  needs  and 
desires  of  the  urban  population  within  their 
ability  to  pay. 

Theoretically,  an  optimal  system  which  best  fits  the 
economic  and  social  objectives  would  be  the  desired  system 
and  the  evaluation  of  the  system  would  be  based  on  criteria 
which  reflect  these  objectives.   However,  the  evaluation 
would  be  very  difficult  if  it  were  to  be  done  quantitatively. 


For  example,  The  Chicago  Area  Transportation  Study  listed 
the  following  six  criteria:  (4) 

1.  Greater  speed. 

2.  Increased  safety. 

3.  Lower  operating  cost. 

It.         Economy  in  new  construction. 

5.  Minimizing  disruption. 

6.  Promoting  better  land  development. 

It  is  apparent  that  a  single  value  system  is  not  available 
for  measuring  all  these  criteria.   Criteria  three  and  four 
are  easily  expressed  in  dollars,  and  one  and  two  can  be 
converted  into  dollar  cost.   however,  criteria  five  and 
six  are  not  only  intangible  but  also  very  difficult  to 
measure  quantitatively.   For  this  reason,  transportation 
system  evaluations  are  generally  restricted  to  economic 
analysis  and  leave  the  social  objectives  to  be  considered, 
somewhat  subjectively,  in  the  final  selection  of  the  op- 
timal system.   It  is,  although  not  perfect,  the  best  ap- 
proach available  at  present. 

This  study  was  an  attempt  to  develop  a  mathematical  model 
for  the  economic  evaluation  of  a  transportation  system. 
Like  many  other  studies,  a  single  objective — to  minimize 
the  transportation  cost — was  adopted.   Several  surveys  have 
shown  that  travel  time  is  dominant  as  a  factor  in  selecting 
a  route  (5,6).   Therefore,  to  simplify  the  formulation,  it 
was  assumed  that  travel  time  was  the  only  significant  factor 


of  operating  cost  thereby  reducing  the  objective  to  min- 
imizing the  suit,  of  the  investment  cost  and  travel  time 
cost . 

One  of  the  major  setbacks  of  linear  programming  type 
models  is  that  unit  travel  time  is  assumed  to  be  independent 
of  traffic  volume.   To  overcome  the  weakness,  this  study 
applied  the  Discrete  Maximum  Principle — a  powerful  optimum 
seeking  method — which  allows  the  use  of  a  non-linear  total 
travel  time  equation.   It  should  be  emphasized  that  although 
there  is  undoubtedly  room  for  improvement  in  the  particular 
functional  relationship  between  travel  time,  assigned  volume 
and  investment  cost,  the  primary  purpose  of  this  study  was 
to  demonstrate  the  usefulness  and  the  ability  of  the  discrete 
maximum  principle  in  solving  this  type  of  non-linear  optimi- 
zation problem.   Hopefully,  by  further  research  and  modifi- 
cations of  the  formulation,  a  more  useful  model  may  be 
evolved. 


REVIEW  OF  LITERATURE 

In  the  past  few  decades,  several  methods  of  economic 
analysis  have  been  developed  for  use  in  transportation 
planning.   Four  principle  methods  are  (a)  annual  cost 
method,  (b)  present  worth  method,  (c)  benefit-cost  ratio 
method,  and  (d)  rate  of  return  method  (7,  &,    9,    10).   The 
relative  advantages  of  each  method  are  briefly  described 
by  Prennan  and  Rothrock  (10). 

Since  World  War  II,  the  benefit-cost  ratio  method  has 
been  given  a  great  deal  of  attention.   In  1952,  the  American 
Association  of  State  Highway  Officials  adopted  this  method 
and  published  an  informational  report  on  "Road  User  Benefit 
Analysis  for  Highway  Improvement,"  the  so-called,  "Red 
Book"  is).      Since  that  time,  it  has  been  accepted  by  many 
planning  agencies.   Grant,  on  the  other  hand,  specially 
favors  the  use  of  the  rate  of  return  method.   The  applica- 
tion of  this  method  and  its  merits  have  been  discussed  at 
length  in  his  papers  (9,  11).   Many  reports  have  compared 
the  uses  of  the  benefit-cost  ratio  method  and  the  rate  of 
return  method  in  detail  (11,  12).   In  general,  when  properly 
used,  both  methods  will  give  the  same  results. 

No  matter  which  method  was  used,  the  analyses  made  in 
the  past  have  restricted  themselves  to  comparing  alternatives 
for  a  single  link  or  a  single  route  of  a  network.   The  over- 
all system  effect  of  improvements  was  completely  ignored  or 
simply  adjusted  by  using  engineering  judgment.   Since  the 


improvements  will  generally  cause  redistribution  of  traffic 
and  since  benefits  on  one  route  may  cause  a  loss  of  benefits 
or.  other  routes,  this  approach  has  sometimes  resulted  in 
an  uneconomic  or  even  retarded  transportation  plan. 

Realizing  this  deficiency,  some  recent  studies  have 
compared  alternatives  through  complete  network  analysis. 
In  the  Chicago  Area  Transportation  Study,  five  alternative 
freeway  systems  were  first  developed.   For  each  system, 
trip  distribution  and  traffic  assignment  techniques  were 
employed  to  obtain  the  traffic  volume  on  each  link  of  the 
system.   Three  methods,  (a)  benefit-cost  ratio,  (b)  rate  of 
return,  ana  (c)  annual  cost,  were  then  used  to  evaluate  each 
system.   They  all  resulted  in  the  same  answer  (4,  13).   This 
approach  is  in  general  quite  satisfactory.   However,  because 
the  number  of  alternatives  to  be  compared  was  relatively 
small  and  the  development  of  alternatives  was  largely  based 
on  engineering  judgment,  it  is  quite  possible  that  the  best 
system  was  not  considered. 

At  the  same  time,  more  and  more  attention  has  been 
concentrated  on  the  applications  of  optimization  techniques 
to  the  transportation  field.   In  1958,  Garrison  and  Marble 
(14)  presented  a  linear  programming  formulation  for  the 
analysis  of  network  improvement.   Travel  cost  for  each  link 
was  assumed  to  be  constant  and  the  investment  was  assumed  to 
increase  the  capacity  linearly.   The  objective  was  to  miminize 
the  sum  of  the  investment  and  travel  cost  subject  to  constraints 


as  flow  balance,  budget,  and  capacity  limits, 
simplex  algoithm  was  employed  to  seek  the  optimal  solution. 
This  paper  leads  to  increased  interest  in  developing  mathe- 
matical models  in  the  following  years. 

Quandt  (15),  in  I960,  presented  a  similar  formulation. 
In  his  problem,  a  commodity  was  to  be  transported  from  n 
sources  to  m  destinations,  each  node  was  connected  with 
every  other  node  by  a  direct  link.   Again,  linear  relations 
between  shipment  cost  and  volume,  and  between  improvement 
cost  and  capacity  were  assumed.   The  problem  was  then  formu- 
lated and  solved  by  linear  programming.   Consideration  was 
also  given  to  problems  with  fixed  budgets,  indirect  con- 
nection between  nodes  and  multicommodity  shipment. 

Carter  and  Stowers  (16),  in  1963,  again  utilized 
linear  programming  to  formulate  a  model  for  fund  alloca- 
tion for  urban  highway  system  capacity  improvement.   The 
basic  formulation  was  the  same  as  the  previous  papers 
except  that  each  link  was  represented  by  two  arcs,  one  with 
free  flow  capacity  and  normal  operating  cost,  the  other 
with  higher  operating  cost  (due  to  congestion)  and  capacity 
equal  to  the  difference  between  possible  and  practical  ca- 
pacity.  The  ratio  of  the  capacities  of  these  two  arcs  was 
kept  constant  as  the  capacities  were  improved. 

In  196/*,  Roberts  and  Funk  (17)  developed  a  linear 
programming  model  for  the  problem  of  adding  links  to  a 
transportation  system.   The  locations  of  possible  additional 


links  in  the  system  were  first  decided.   In  seeking  the 
optimum,  the  additional  link  was  either  completely  built 
or  not  built  at  all.   If  the  link  were   ided,  the  cost 
was  included  in  the  objective  function.   If  it  were  not 
added,  the  flow  was  blocked.   In  this  formulation  an 
integer  programming  technique  was  used.   The  paper  also 
suggested  a  possible  application  of  dynamic  programming 
in  treating  the  stage-wise  construction  problem.   As  a 
result,  in  1966,  Roberts,  et .  al.  (18)  combined  the  use 
of  linear  programming  and  dynamic  programming  techniques 
to  solve  a  stage-wise  link  adding  problem.   The  budget  for 
each  construction  period  was  fixed.   The  method  considered 
the  budget  at  the  Nth  stage  as  the  sum  of  the  budgets  from 
the  first  to  the  nth  construction  period  and  used  the  prin- 
ciple, "the  links  for  stage  K  must  be  the  subset  of  links 
for  stage  K  +  l,1'  to  indicate  the  links  which  must  be  con- 
sidered at  stage  N.   At  each  stage,  integer  programming 
was  used  to  select  the  links  to  be  added.   This  method  is 
considered  useful  for  transportation  system  development  in 
underdeveloped  countries . 

In  the  same  year,  Kay,  Korlok  and  Charnes  (19)  presented 
a  model  for  optimal  planning  of  a  two-mode  urban  transpor- 
tation system.   A  two-mode  system,  private  road  and  public 
transit,  was  to  be  built  in  an  urban  corridor.   The  road 
capital  cost  was  linearly  related  to  capacity  and  speed. 
Transit  speed  was  fixed  with  the  capacity  linearly  related 


to  capital  cost.   The  length  of  the  transit  route  was  also 
assumed  fixed.   The  choice  of  mode  was  linearly  related  to 
the  travel  time  ratio  between  road  and  transit.   Again,  the 
linear  programming  technique  was  used  to  formulate  the 
problem  and  seek  the  optimum.   The  objective  function  was: 

fannual  road  ]   fannual  vehicle]  ("annual  transit 

mirM  H  M    ,   , 

(capital  cost!   (operating  costj  [capital  cost 


annual  transit 


\operating  costj +<Parking  costj 


In  this  formulation,  the  travel  time  was  excluded  from 
operating  cost  and  was  treated  as  a  constraint  to  reflect 
the  minimum  level  of  service  desired  and  the  maximum  speed 
obtainable.   For  a  true  optimum,  it  was  required  to  change 
the  length  of  the  transit  route  and  run  the  program  several 
times . 

Distinct  from  the  linear  programming  type  models, 
Ridley  (20)  in  1965,  developed  a  method  for  seeking  the 
optimum  investment  policy  to  reduce  the  travel  time  in  a 
transportation  network.   The  unit  travel  time  was  assumed 
to  be  decreasing  linearly  with  the  investment.   It  was  also 
assumed  that  the  flow  was  far  below  the  link  capacity.   The 
objective  was  to  minimize  the  total  travel  time.   Because 
the  travel  time  was  a  function  of  both  investment  and  traffic 
volume,  the  objective  function  was  non-linear  in  nature.   For 
some  special  cases,  such  as  no  budget  limit,  fixed  traffic 
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volume,  fixed  investment  and  single  origin-destination, 
the  formulation  can  be  simplified  into  a  linear  program- 
ming model.   For  the  general  case  which  has  budget  and 
travel  time  constraints,  the  bounded  subset  method  was 
utilized  to  search  the  optimum. 

A  common  drawback  of  the  above  models  is  that  unit 
travel  time  was  assumed  to  be  independent  of  traffic  volume 
which  is  not  at  all  close  to  reality.   The  linear  assump- 
tion of  travel  time  improvement  in  Hay's  and  Ridley's 
models  is  also  questionable. 

Although  the  above  review  shows  some  imperfection  of 
the  existing  models,  it  is  realized  that  the  mathematical 
model  is  a  powerful  technique  which  we  can  apply  in  the 
field  of  transportation  engineering  in  order  to  make  a 
better  analysis.   Compared  to  the  alternative  comparison 
method,  the  mathematical  model  has  merit  in  that:   (19) 
''It  selects  that  system  which  is  optimal  among  all 
possible  systems  of  a  given  type  rather  than  merely 
examining  a  small  number  of  alternatives.'' 
There  are,  however,  two  difficulties  with  the  mathematical 
models. 

1.  We  are  trying  to  express  a  highly  complex  real 
world  phenomenon  with  a  simple  equation  which  is 
very  difficult  or  may  even  be  impossible. 

2.  The  model  sometimes  becomes  too  complex  to  man- 
ipulate and  too  sophisticated  to  understand  and 
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requires  highly  trained  personnel  to  carry  cut 
the  analysis. 
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;te  maximum  principle  and  its  recent  APPLICATIONS  IK 

TRANSPORTATION  SYSTEMS  ANALYSIS 

In  recent  years,  optimization  has  become  more  and  more 
important  in  engineering  systems  analysis-   Since  World 
War  II,  many  sophisticated  techniques  have  been  developed. 
Among  them,  Dynamic  Programming  and  the  Discrete  Maximum 
Principle  (1)  were  developed  for  the  optimization  of  stage- 
wise  processes. 

In  196/,,  ran  and  Wang  developed  a  discrete  version  of 
the  maximum  principle  (1).   Recently,  it  was  demonstrated 
that  this  optimization  technique  was  applicable  to  tran- 
sportation systems  analysis  (22,  23,  24). 

Consider  an  N  stage  process  with  state  variables  de- 
noted by  an  s-dimensional  vector,  X  «  (X  ,  X?)...  X  ),  and 
decision  variables  denoted  by  a  r-dimensional  vector,  9  = 
(a-,  ,  89,... Sr).   The  performance  equations  at  the  n-th  stage 
are  given  as: 

x»-T»(iri.xr1.---.*r1i  *?.•§.-.# 

where,  i  =  1,2, ...,s;  n  =  1,2,. ..,N  and  &j_  is  constant. 

A  typical  optimization  problem  associated  with  such  a 
process  is  to  find  a  set  of  6n,  n  =  1,  2,...,  N,  subject  to 
constraints : 
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ty  j  (0y,c?,...,3")  —  °         n  =  1,  2,...,  K 

i  =  1,  2 r 

which  optimizes  the  ohjcetive  function, 

s        ., 
S  =  2Z  Cj  x  X.  C.  =  constant 

The  discrete  maximum  principle  introduces  an  s-dimen- 
sional  vector,  Zn,  and  a  Harniltonian  function,  H  ,  which 
satisfy  the  following  relations: 

Kn  =  2Z  Z?  Tn  (Xn_1;  on),     N  =  1,  2 N 

i-1  *   1 


Z""1  =  -2 — r      i  -  1,  2,...,  s;   n  -  1,  2,...,  N 

xn-l 

w 

and  Z .'  =  C  . 


The  necessary  condition  for  S  to  be  a  local  extreme  with 
respect  to  S  is 
°H =  0 

when  it  is  inside  the  boundary  of  the  constraints,  or 

Hn  =  extreme 
when  it  is  on  the  boundary  of  the  constraints. 

If  the  objective  function  is  a  cumulated  measure  which 
can  be  expressed  as 


n-1  ' 
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the  al  :o)  i  1  a   car.  be  extended  by  introducing  an  extra  state 
iefined  as 

XS+1  "  ° 

vn     ,,n-l   ,|.  ,,,n-l   n, 

xs-rl  =  xsil  -1  Y  (X   ;  a  ),     n  -  1,  2,...,  N 

new  equations  together  with  the  original  performance 
equations  specify  the  process  in  s+1  variables  where  the 
objective  function  becomes 

n=l 

and  the  primal  algorithm  of  the  principle  is  restored. 

In  case  that  some  of  the  state  variables  at  the  end 

N  N 

stage,  Xi,  are  fixed,  the  relation  Z±  =   Ci  no  longer  exists. 

The  following  equation 


can  be  used  for  each  fixed  end  stage  variable.   By  solving 

N 

these  equations  simultaneously,  the  Z.  values  can  be  determined 

(21). 

Recently,  the  discrete  maximum  principle  has  beer,  suc- 
cessfully applied  to  the  traffic  assignment  problem.   In 
1964,  Yang  and  Snell  (22)  formulated  the  traffic  assignment 
problem  by  considering  each  node  of  the  network  as  a  stage. 
The  unit  travel  time  on  each  link  was  assumed  constant  and 
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a  turning  oer.alty  was  included  in  the  travel  time  equatior. . 
The  objective  was  to  minimize  the  total  travel  time.   Due  to 
the  linear  characteristic  of  the  objective  function,  the 
optimum  searching  procedure  was  reduced  to  a  shortest  path 
tree  building  routine  which  is  similar  to  Moore' 6  algorithm. 

In  1966,  Snell,  Funk  and  Blackburn  (23)  developed  a 
more  complete  model.   In  this  model,  travel  time  was  non- 
linear ly  related  to  traffic  volume.   This  characteristic  is 
considered  to  be  a  step  toward  more  realistic  traffic 
assignment.   In  1967,  Funk  and  Snell  (2/,)  developed  a 
procedure  for  an  approximate  multicopy  traffic  assignment 
problem.   The  results  obtained  from  this  procedure  appear 
to  be  very  close  to  the  true  optimum. 
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THE  objective  function  and  travel  tike  equation 

As  previously  explained,  the  objective  of  this  study- 
was  to  minimize  the  sum  of  the  investment  cost  and  the 
travel  time  cost.   The  investment  was  an  independent  var- 
iable and  it  was  assumed  ihat  it  could  be  expressed  in 

.'  of  dollars  per  mile.   However  unit  travel  time  was, 
Ln  ;eneral,  dependent  on  traffic  volume  and  roadway  con- 
ditions.  In  other  words,  unit  travel  time  was  a  function 
of  both  traffic  volume  and  investment.   The  relationship 
among  them  was,  in  reality,  very  complex.   In  developing 
a  mathematical  model,  it  was  generally  necessary  to  make 
some  assumptions  and  simplify  the  relationship  in  order 
to  express  the  relationship  by  a  relatively  simple  equation 
which,  was  manageable  and  yet  not  too  far  from  reality. 

To  express  unit  travel  time  as  a  function  of  traffic 
volume  and  investment,  some  basic  characteristics  were 
observed : 

1.  Unit  travel  time  increased  as  the  traffic  volume 
increased. 

2.  Unit  travel  time  decreased  as  the  investment  in- 
creased . 

3.  Unit  travel  time  had  a  lower  limit.   (free  flow 
travel  time) 

4.  With  constant  travel  time,  service  volume  increased 
as  the  investment  increased. 
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The  typical  relation  between  traffic  volume  and  oper- 
ating speed  is  shown  in  Fig.  la.   As  the  speed  is  inversely 
proportional  to  the  travel  time,  this  curve  can  be  converted 
into  a  travel  time-traffic  volume  relation  curve  as  shown 
in  Fig.  lb.   The  dotted  part  of  the  curve  shows  the  relation 
under  congested  conditions.   Therefore,  under  normal  oper- 
ating conditions,  it  is  logical  to  assume  that  unit  travel 
time  (in  hours  per  vehicle  per  mile)  is  linearly  related  to 
traffic  volume  and  should  have  an  equation  of  the  following 
form : 

t  =  K  +  K'V  (1) 

where 

t  =  unit  travel  time  (hr/mi/veh) 

K  =  free  flow  travel  time  (hr/mi/veh) 
K'  =  slope  of  the  curve  in  Fig.  lb  (hr2/mi/veh2 ) 

V  =  traffic  volume  per  unit  time  (veh/hr) 

Keeping  basic  characteristics  in  mind  and  further  as- 
suming that  the  free  flow  travel  time  is  constant  for  each 
link  and  traffic  volume  served  is  proportional  to  investment 
for  a  constant  travel  time,  an  equation  of  the  following 
form  may  be  hypothesized: 

t   =  Kx  +  hv  (2) 

where 

t  =  unit  travel  time  (hr/mi/veh) 
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Link  Volume ,  vph 
Fig.  la  Typical  Speed-Volume  Curve 


Link  Volume,  vph 
Fig.  lb  Typical  Travel  Time- Volume  Curve 
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K,  =  free  flow  travel  time  (hr/mi/veh).   The  magnitude 
depends  on  the  maximum  speed  obtainable  or  regu- 
lated. 
Kp  =  coefficient  of  improvement  (dollar-hr/mi  /veh  ). 

Its  magnitude  depends  on  link  location  and  reflects 
the  difficulty  of  improvement. 
8  =  equivalent  hourly  investment  per  unit  length 

( dollar/mi /hr). 
V  =  traffic  volume  per  unit  time  (veh/hr). 
In  the  case  where  old  facilities  exist,  the  investment 
should  be  expressed  as: 

(3) 

where,  Kj  ,  in  dollars  per  mile  per  hour,  represents  the 
existing  investment  and  91 ,  in  dollars  per  mile  per  hour, 
is  the  additional  investment. 

The  general  form  of  the  unit  travel  tine  equation  then 
becomes 

t  =  jC  +  __! r  V  (4) 

The  characteristics  of  this  equation  are  demonstrated  in 
Figure  2a,  2b  and  2c. 

Let  L  be  the  length  of  the  link  and  Ct  the  cost  of  time. 

The  objective  function  then  becomes 

S  "  9'L+(K1V+  ±1^   V2)  L  Ct  (5) 
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b<b„<i, 
12  3 


Investment,  6 
Fig.  2a  Travel  Time-Invest- 
ment Curve  With  Fixed 
Volume 


Traffic  Volume,  V 
Fig.  2b  Travel  Time-Volume 
Curve  With  Fixed 
Investment 
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Traffic  Volume,  V 
Fig.  2c   Volume-Investment  Curve 
With  Fixed  Travel  Time 


21 


GENERAL  FORMULATION  OF  THE  PKORLEM 

This  section  presents  the  general  formulation  of  the 
optimal  network  improvement  problem.   Three  investment  con- 
ditions are  considered  which  result  in  three  different 
sets  of  equations  and  two  slightly  different  methods  of 
seeking  the  optimum. 

The  formulation  is  aoplied  to  rectangular  networks  only. 
However,  for  many  non-rectangular  networks,  it  is  possible 
to  modify  them  into  rectangular  forms  by  adding  slack  links. 
Two  examples  are  shown  in  Fig.  3. 

Definitions 

1.  Objective  Function  -  a  function,  which  is  to  be  minimized 

in  this  problem,  representative  of 
the  total  cost. 

2.  Zone  Centroid  -  a  point  of  trip  origin  or  destination. 

3.  Node  -  a  point  where  segments  of  the  road  system  connect, 
i .   Link  -  a  connection  between  two  nodes  representative 

of  a  segment  of  the  road  system. 

5.  Path  -  a  series  of  connected  links  representative  of  a 

trip  route. 

6.  Network  -  the  combination  of  all  links  and  nodes. 

Notations 

1.   X.'   -  state  variables  representing  flows  from  node 
(n,m) . 


^> 


=> 


-(•>- 


-(■)- 


slack  links  with 
zero  travel  time 


slack  links  with 
infinite  travel 
time 


Fig.  3  Examples  of  Network  Modification 


n  in 

2.  6j'   -  decision  variables  representing  investments  on 

links  leaving  node  (n,m). 

3.  Kj£  -  Free  flow  travel  time  on  links  leaving  node 

( n ,  m ) . 

4.  K^Q  -  coefficient  of  investment  on  links  leaving  node 

(n,m). 

5-   K-4   -  existing  investment  on  links  leaving  node  (n,ra). 

6.  LV'1*  -  link  length  on  links  leaving  node  (n,m)  . 

n  m 

7.  tj   -  unit  travel  time  on  links  leaving  node  (n,ra). 

where,  j  =  1,  for  horizontal  links. 
j  =  2,  for  vertical  links. 

v    i-n  m 

5.  X,  "  -  state  variable  representing  the  total  investment 

on  horizontal  links  from  node  (1,1)  through  node 
( n ,  m ) . 

9.   X, '   -  state  variable  representing  the  total  investment 
on  vertical  links  from  node  (1,1)  through  node 
( n ,  m ) . 

10.   Xr' '  -  state  variable  representing  the  total  travel 
time  cost  on  horizontal  links  from  node  (1,1) 
through  node  (n,m). 

11-   x6''  "  5tate  variable  representing  the  total  travel 
time  cost  on  vertical  links  from  node  (1,1) 
through  node  (n,m). 


.  I, 


12.  xS,m  -  state  variable  representing  the  total  investment 

on  both  links  from  node  (1,1)  through  node  (n,m), 

13.  d!?'^  -  decision  variable  representing  the  fraction  of 

the  vehicles  departing  node  (n,m)  on  the  hori- 
zontal link. 
Lrr.e  cost. 

15.  Kn,n  -  Hamiltonian  function  at  node  (n,m). 

16.  V  >m  -  input  trips  at  node  (n,m). 

17.  GI  -  total  system  budget. 

IS.   SIn,r'  -  section  budget  at  node  (n,m). 

19.  Z?,m,  Z^'"1,...,  Z-"',r-  -  adjoint  variables  associated 

with  Xy»  ,  X2   ,.•-,  aj 
respectively. 

20.  S  -  objective  function. 

Assumptions 


1.  Ko  turn  penalties. 

2.  Zone  centroids  coincide  with  the  nodes. 

3.  Traffic  directions  are  preassigned. 

4.  Traffic  distribution  is  fixed. 

5.  Transportation  network  can  be  represented  by  a  rectangu- 
larly arranged  combination  of  links. 

6.  Travel  time  is  the  only  factor  that  influences  the 
traffic  assignment. 


2  J 


7.    Unit  travel  time  on  each  link  can  be  expressed  as: 


"J      Jl    en'm  H  K^"1"'1 


Kn,m 
12 xn,m  (6) 


where 

j  =  1,  for  horizontal  links 

j  =  2,  for  vertical  links 

Figure  U   shows  a  basic  N  x  M  rectangular  network  with 
node  (N,M)  as  the  destination  and  all  other  nodes  as  origins. 
With  the  input  trips  at  each  node  obtained  from  a  traffic 
distribution  study,  the  problem  is  to  find  an  investment 
policy  under  each  investment  condition  such  that  the  total 
cost  is  a  minimum. 

Investment  With  No  Budget  Constraint 

In  this  case,  the  overall  system  budget  was  assumed 
unlimited.   However,  there  are  three  special  conditions  which 
imply  upper  or  lower  limits  of  investment  on  each  link. 

The  performance  equations  for  a  typical  interior  node 
as  shown  in  Fig.  5  were  developed  as  follows: 

xn,m  ,  (xn,m-l  H  xn-l,m  +  vn,m)en,m  .  AIn,men,m  (?) 

xn,m  =  (xn,m-l  +  xn-l,m  +  ^.m,  (1.eM,  .  Al^d-a"*")   (8) 
xn,m  =  eri,m  Ln,m  +  xn,m-l   ;    en,m  >0  (9) 

x"-n  =  e"'m  L5'm  4  x"'"1-1   ,    e"'m  >o  (io) 


26 


VK-1,M- 


TnTD  (ImT 


1*1^  * 


i  t 


(N-l.M) 


(N.M-1)  \J[N,M) 


Fig.    4      Basic   N  x  M  Network 
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n-1 


n+1 


m-1 


(n-l,m-l) 


Yn,ia-1 


(n,m-l) 


(n-l,m-I) 


(n-l,m) 


„n-l,m 
A2 


rn,m 


m+1 


(n,m) 


Yn,m 
X2 


(n+l,m) 


(n-1, m+1) 


( n , m+1 ) 


(n-rl.m-rl) 


Fig.  5  Typical  Interior  Node  of  A  Rectangular 
Network 


28 


,  n ,  m   r  n ,  rn 


vn,m        -n.ir,   vn,ra   Tn,m  r       ,      12       1  t    /Yn>nM2  +yn,m-l 

K5        '    Kll       '1  1  t    '    ^,n:    ,    Kn,ra       Ul       '     'X5 


B1         I    K13 


„n,m   ATn,m„n,raTn,r.„  "'12     Jjl        ut    .      n,m  n,ni,2 

°1        "*    I'll 


13 
n,m-l  (ID 


„n,m  _  „a,B  yn,m  Ln,m  .     ,  K22 '    V        ct       n,m ,2       rn,m-l 

X6       "  -21     A2        L2       Ut     cn,ir.       rii,n        (x2      '      ^x6 


rn,m       n1a/1_nn,ii,Tn,!iiw 

"  +  K 


n,n  Tn,m 

22 

n, 

2 

vn,Ei  Tn,m  „ 

n,n 
23 


M    an,m,Tn,mr      ,    ^22     ^2       nr.n,a,      hn>mj|2 

[i-e3«  )l2>  ct  +  j  vJiiU  |ai      (i-e3    j 


_   yn.m-l 

-6'  (12) 


where  ^n.m  =  :(n,r.-l  +  xn-l,m  +  ^.n  (13  } 


=2  0 
>0 


ind      0  <  Sj'"^  1 

The  Kamiltonian  function  at   this   node   is   defined  as: 


-n,^_    7n,m  vn,m        „n,m  „n,in    ,    „n,r.  vn,rn    ,    „n,m  Yn,n 
-  i-j_        A-]_        +    l2        x2         i-  ^3        Xj         +   ^        A^ 


+  2».»  X^  +   2g.»  Xg."  (14) 
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.     Ltutinp   equations    (7)    to    (12)    into   equation    (1/J 
and    takir./i   derivatives   with   respect   to   state   variables,    the 
adjoint   variables  are   obtained  as   follows: 


n.m-1  _    3  Hnim       „n.nun.m    ,    „n,m/,    rtn,m> 

1  "a^T^T  "     1     b3       +  Z2      (1-a3      ' 


8X 


,    „n  ,iuvn,ni,  n,mTn,m- 
+   Z,      K, ,'    9-  *    L-,  '    C 


'1 


7n,mKn,m;1    w   n,rr.  i    Tn,mr 
Z6     K21    (1-b   3>    )    L2'    Ct 


n,m  vn'n  t"'131 


V      +  Klj 


.  _n,m,.,n,m,  2 
AI        (b3'    ) 


v22 


n,m/,_ fin,m  2 


,n,m       vn,in 


t    Ain'm(i-e3"Ui) 


v23 


(15) 


r.n-l,ir,  _    3 


V  =  ?n,mfin,iri  +  7n)rifi    «n>mi  4.  7n,i%n,mQn,mTn, 

r.-l,rn       ll     b3  V    (1_ri3      '   T  z5     Kli  d3     Ll 


6      ]121    U   d3       '      2        t 


n,m        n,mr 

h   2  z"'m  "12     Li h    ai— (8?"")' 

2       an,m    ,    B-n.m  3 


-n.ffl,   n,m.2 


>n,m 

31 
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Kn,m   Ln,mc 
zn,n;  _22 2 1  A1n,m(1_en,ni)2 

6        „n,m  n,rr.  A1        u   °3      j 

b2        ~  K23 


(16) 


7n,ri-±  =dh    '  '  r7n,m 

"3  "d^^    "     '3 


(17) 


30 


7n-l,m  =  3  i:r''m  =  zn,m  (lf>) 


-  dxn-l,m 

4 


'5     ~  I^1   "  ^ 


zn-l,m  „  5  Hn'"'  =  zn,m  (20) 

6       iyr.-l,ra    6 


The  original  conditions  for  the  state  variables  are  given  as: 

xo,o  =  xo,o  m   xo,o  m   xo,o  =  xo,o  =  xo,o  =  Q  (21, 

The  objective  function  is 

„N,M   „N,M   VN.M   „N,K  ,„„» 

S  =  X3-   +  X,'   -  X5'   +  X6«  (22) 

Therefore,  by  definition,  the  boundary  conditions  for  the 
adioint  variables  are: 


z^f1-  zV"  -  0  (23) 

zN,v=  ZN,K=  ZK,M=  ZK,K=  x  (24) 


Substituting  equation  (24)  into  equations  (17)  to  (20),  the 
following  ecuation  is  derived. 


3-1 


.  ::f'n!  -   Z"'nl  -   Z?'m  -  1         for  all  (n,m)     (25) 


The  Hamiltonian  function  then  becomes: 

Hn,m  .  zn,mxn,r.  4  ^.n^n.m  +  xn,m  +  ^n.m  4  ^.m  +  xn,m   (2fi) 

Ir.  order  to  have  S  a  minimum,  the  following  conditions  are 
necessary : 


IK „  0  en,m  >0 


dH^f  =  0  en,m  >0 

a8n,r.  2 


^f  =  0  o<en>m<i 

se">ra  3 

/  *i  jpd   r1  m   r  n  ^  rim 

when  ISi'  ,  ct"'  ,  bJ'  '  J  is  an  interior  point,  or  H    ■ 

n  ,m 
minimum  with  respect  tc  those  «■   which  are  at  a  boundary- 
point  of  the  constraints. 

Substituting  equations  (7)  to  (12)  into  equation  (26) 
and  taking  derivatives  with  respect  to  the  various  decision 
variables,  the  following  equations  are  obtained: 


,..r.,m  K"'m  L?>m  C 

olj __  Tn,m    !<;   1    t      n  m  n  m  2 

..n,m  ~  ^1    "   n,rr.     ,n,m  2  (AI  '  ei'     )  (27) 


3  „*!•*«         -L  .       ll.il.  li  .  lil 

°1  («i    T  K13  ) 
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dh ,n,m  ,:.:  r.Tn,m,        n.m,-.^  ,„,.> 

^Ti        L2        "    (onlK   "    Knf)2   [XI         (I-B3       )}  (26) 


/7n,iT.        7n.m>.Tn,Bi        /-.n,mTn,m  v.n,mTn,m> >Tn,nu 
(Zx>      -   Z2>    )A1  -    U>    L]_>    -K2i   L/    )AI        Ct 


00,  ■ 


n ,  ra     n ,  rr. 

2Hr.,.T.  ^   Kr.,ir.        (A1         '    °3 
«  «   K13 


r*  m     r  rr* 

Taking  derivative  of  equation  (29)  with  respect  to  8  '  ',  the 
following  equation  is  obtained: 


_  2_12 1 1  {AIn,m}/  +  2_J2 2 1  (AIn,m)2  (30) 

,n,n.  ,  ,.n,m  n.m   vr.,m 

dl   4  K13  °2   +  K23 


Setting  equations  (27)  and  (28)  equal  to  zero  and  apply- 
ir.r  '-he  boundary  conditions  of  the  decision  variables,  the 
values  of  d ''   and  9  '  ''  can  be  obtained  from  the  following 
eauations : 


ei'm*yK12mCt  AI0'"^'"1  -   K^m     when  d^>rr->0  (3D 
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n  .:"■  /,  n,ir,„  ■"   ..n,^  <--  „  i„  i 

hi'  -   U    when  ,/K12  Ct  A,.   83   -  KjJ   <0  (32) 


8n,m=  /K^lr-Ct  AIn'm(l-en'm)-  Kr2f     when  8°.»  >0        (33) 


6,  =  U  wnen  V^22  t      (I-O3   J-  K-23  —  l>^' 

When  both  e",m  and  e2>m  are  greater  than  zero,  equations 

(31)  and  (33)  can  be  substituted  into  equation  (29)  to  obtain 
the  following  equation 


_  ,  7n,m  7n,m.  ATn,m  _  ,   n,m  n,m  n,m  n,n,   n,r.c 
„  n  a  *=  I  zl  "z2   '  AI     •  (  Kll  ^1   K21  ~2   'Ai   ut 


„  Opu"  Tn,m,Tn,m   „  OTTnZ"  Tn,m.Tn,Ei 
n  2/K^  C^'  AI    -  2  ,/K22  CtL2  Al 

.  AIn,m  ((  zn,m_zn,m)  +  (  ^^n.m^mq,^ 
+  2t/iff5;  I-.-  .  /ign^  i».»,)  (35) 


.n.m 


60'   is  eliminated  by  the  substitution  and  the  value  oi 

n,n; 

^ "n  ^  becomes  independent  of  83 '  '  as  shown  in  equation 

n  ni 
(35).   This  ir.piies  that  the  value  of  H    is  linearly 

related  to  9,   ana  the  extreme  of  H    with  respect  to 

83 '  ':  occurs  at  a  boundary.   In  this  case,  to  obtain  the 

,n,m 


,_•'-.- 


i 


3/. 


When  I^t— n:  <^   If  '  r'i-'  =  °>  D3,K1  can  be  any  value  between 

n  m  n  > lr* 

0  and  1  because  the  value  of  H  '   is  independent  o^  Bj 

When  either  8?,m  or  o^'m   is  equal  to  zero,  or  when 
both  are  equal  to  zero,  equation  (35)  is  no  longer  valid. 
Equation  (29)  is  then  set  equal  to  zero  and  solved  for  the 
optimal  value  of  6,   . 

Special  Case  I:   In  an  urban  area,  the  available  space 
for  road  construction  is  often  limited.   For  example,  a  free- 
way with  more  than  eight  lanes  would  be  very  difficult  to 
build  near  a  CBD  area.   It  is,  therefore,  necessary  to  set 
an  unper  limit  on  the  size  of  the  links.   These  limits  can 
be  expressed  as  limits  on  investment.   Mathematically,  they 
are  expressed  as: 

,,n,n  ,  ,n,m  <--  „n  ,m 
Kll   "■  bl   -°1  max. 


R21  n  b2   -92  max. 

Equations  (3D  and  (33)  are  then  replaced  by: 


,-,-n.m   ^n,m   ,,n,m 
d3   "  A13 


^1'max.  "  Klim»  (36) 


>n,m    /„n,nu   ,  Tn,mn  ,n,n,   „n,r 


The  other  equations  remain  unchanged. 

Special  Case  II:   In  developing  a  urban  transportation 
network,  it  is  sonetir.es  required  to  provide  a  minimum  level 
of  service  for  the  entire  area.   For  example,  arterial  streets 
would  be  distributed  uniformly  throughout  the  whole  area. 
This  criterion  can  be  fulfilled  by  requiring  a  minimum  amount 
of  investment  on  each  link.  Mathematically,  it  can  ,oe  ex~ 
pressed  as: 

Kn,m  +  en,m  ^n.a^  (38) 

,.n,m  j  „n,m  ->an,m  (to) 

K2i    °2        — d2  min. 

Special  Case  III:   When  the  conditions  governing  both 
special  cases  I  and  II  exist,  both  upper  and  lower  limits 
of  investment  should  be  applied  to  each  link.   Mathematically, 
they  are  expressed  as: 

Qn,m    <-  „n,n   ,.n,m<,,n,m  /,nl 

8lmin.-Kli  T  °1   -°lmax.  (40) 


rr.,m   <-  „n,m  ,  Qn,m  <- Qn,m  /,,i 

e2  min.-  K2i  4  V  -  b2  max.  (U) 
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•  above  formulation  provides  the  equations  for  search- 
ing the  optimum  sequence  of  the  decision  variables,  6^'w', 

gn,;r,  an(j  ^5J>m  ariCj  the  associated  values  of  the  state  var- 
2        J 

iables . 

The  optimum  seeking  procedure  developed  for  this  problem 
is  as  follows: 

1.  Assume  a  set  of  decision  variables,  8_j'  . 

2.  Calculate  X",m,  X?  '  "  and  AI  '   by  equations  (7) 

1     *- 

(8)  and  (13)  starting  at  n=m=l  and  proceeding  to 
n=N  and  m=M. 

3.  Calculate  decision  variables,  o-,'  '  and  9,'  ,  by- 
equations  (31)  and  (33)  and  check  the  boundary 
conditions  for  each  special  case. 

k.  Calculate  the  values  of  X^'m,  i  =  3,  4,  5,  6,  by 
equations  (9)  to  (13)  starting  at  n=m=l  and  pro- 
ceeding to  n=N  and  m=K, 

n  ,m 

5.  Calculate  the  adjoint  vectors,  Z_.   ,  i  =  1,  2,  with 

the  above  X- '  "   values,  by  equations  (15),  (16), 
(23)  and  (25)  starting  at  n=K,  m=M  and  proceeding 
backward  to  n=m=l. 

6.  Using  the  above  values  of  X. '  '  and  Z?,m,  calculate 

dHn,m      a2i^-,rr. 

^ ">.,^  and  , .  n   n  2     by  equations  (29)  and  (30). 

7.  Adjust  the  values  of  0?,m  by  adding  an  amount  equal 
to  A ,  where 
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A. 


3;— 


■\2„n  ,m 

o  ■- 

n  m  2 
06,'") 


^3 

and  check  the  boundary  condition. 

n  m 
S.   With  the  new  values  of  8>   ,  return  to  step  2  and 

repeat  the  procedure  until  the  value  of  the  ob- 
jective function,  equation  (22),  is  sufficiently 
close  to  the  previous  value  to  indicate  adequate 
convergence. 

Investment  With  Fixed  Node  Investment 

In  developing  a  large  area  trunk  line  system,  the 
budget  for  each  traffic  section  is  sometimes  predetermined. 
Suppose  that  we  consider  each  node  and  its  associated  two 
links  as  a  traffic  section  where  a  fixed  budget  is  allocated, 
then  a  different  formulation  could  be  developed. 

The  budget  condition  in  this  case  can  be  expressed  as: 

en,m  +  en,m=  £:n,m  a,} 

where,  o't1'1"  and  9?'   are  total  investments  instead  of  invest- 

n,m 
ments  per  unit  length  as  in  the  previous  case.   Let  B\       be 

n  m 
the  independent  variable,  then  d?   can  be  expressed  as: 

6r.,m  =  SIn,m  _  Gn,m  {^} 
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Since  the  total  investment  is  fixed,  the  individual 
investment  costs  are  no  longer  included  in  the  objective 
function  and  also  will  not  be  expressed  as  state  variables. 
The  performance  equations  for  a  typical  interior  node  as 
shown  in  Fig.  5  can  be  written  as  follows: 

xn,m  _  (xn,m-l  +  xn-l,ra  +  ^.m^m  _  ^.m^.a         (u) 
X*'m  -    (X^-'"1  +  X^1*18  +  Vn'm)(l-9^m)  =  Aln'm(l-b5'ir')  (45) 


„n,m  TnJra  p 
vn,m   „n,m  Tn,m„    .Tn,m  an,m  _^  12   1    t  ,  .Tn,mflr.,m,2 
V  =  Kli  V  ct  Ai    e3   T  ^TS        UI   e3  ' 

Tn,m  +  K13 

+  X"'"1"1  (46) 


5 


-■6   =  K2l  Ll  ct  AI   (1-^3   ) 


[Ain'B(i-ef'"))2+'.^;B-1    (47) 


t  r;  i m  r 


,Tn,m  ,,n,ir. 


rn,m     23 


n,r."i       „n,m-l       ,,n-l,ni       ,,n,m  i,r>\ 

where,    AI    *      =  X    '  -.    X2      '      +  V    '  (43) 


0<er.,m  <x 


ana     0    <  e£'m<  SIn'E 


39 


.  amiltonian  function  in  turn  becomes 

n,m   ,r.,:n  n,in  .  7n,ra  Yn,m   „n,m  Yn,m   „n,m  Yn,m  ,   . 
H    -  ^   >-x   h  ^2   X2   h  Zj   Xj   i  Z6   Xg   (49, 

the  cb.lective  function  is 

S  =  X^'M  +  Xp:  (50) 

The  values  of  the  adjoint  variables  are  as  follows: 
zn,m-l  =  zn-l,a  .  zn,men,m  +  zM(1.eM)  +  K^mL^r'Ct  (l-e£'m) 


j(  n ,  in  ^n ,  m  q 

+  2  -12 1 1  AIn,m(en,mj.a 

9n,iB  +  Kn  m        3 

"rTTm 
Ll 


+  2  ;2,  2   t    Ain'B(i^en'm)2         (5D 

SI   -8-,   -  K2,  •? 


2n,m  =  zn,n  =  x  for  &11  (n>m)  ($2) 

and 

Z1'  =  Z2'  =0  (53) 

The  initial  values  of  state  variables  are 

>:o,o  .  xo,o  .  Yo,o  .  Yo,o  =  Q  (54) 


^0 


To  find  the  minimum  value  of  S,  the  following  conditions  are 
necessary 

d-¥^  =   G  0  <6?>m<SIn'r;i  (55) 


^S=0  0<b^r'<l  (56) 


when  I  6,*'  '  ,  6  "'   is  inside  the  boundary,  or 


'X   '  -3 


(57) 


when   [S^'  ,  83'   lis  at  a  boundary  point  of  the  constraints. 

Substituting  equations  (kU)   to  (48)  and  (52)  into 

equation  (A9)  and  then  taking  derivatives  with  respect  to 

9  '"'  and  6n,:K,  the  following  equations  are  obtained: 
l       3 

,  n,m        Kn,a  c 

3K_  m   _   ^-2   t Uln.men>m)2 

ael'      ,_V_    n,m)2 
'T^Tm  +  K13  ] 

"1 


c  n,mK2n,mCt ~     [Al^d-b^) }  *  (58) 


1     _r.,m   T  K23  J 


u 


■,,,n,rn        v   i  2  11       1  <2l       <d         t 

5 


Kn  re   Ln,m  c 


l< 


n,m 

£l T   Kn,m 

rii, in         13 

"1 


t    Urn,m^      fln,m 


(Al"'m) 


n,rr.     n  ,m 
K22     L2        ct 


(AIn'm)    2    (l-e1?'111) 


n,m 

l2 


Kn,m 

23 


(59) 


Setting  equation  (58)  equal  to  zero  and  solving  for  y 
we  obtain 


n,m 


/,.n,m,„Tn,iii  T,n,m,n,m.    n,mTn,m  Ln,m..n,m ,.    ,.n,m,  Tn,mTn,m 
3n,r:_^12    (SI        "k23    L2      >e3      Ll      -jft22_Ki3    (1~d3      )L1      L2 
1  /  ,.n,rr. , ,      ,r.  ,m  |  T  n,m        Xn.n   ,,r.  ,rn    rn,rr. 

V    K22    '1_fcl3       >L2  ^K12      D3  1 


(60) 


Using  equations  (UL)    through  (53),  the  optimum  seeking 
procedure  developed  in  the  previous  case  was  again  applied 
to  solve  the  Droblem. 
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Investment  V.  :t.l-  Fixed  System  Budget : 

It  is  not  unusual  for  the  total  budget  for  a  transpor- 
tation system  improvement  to  be  predetermined.   In  this  case, 
the  total  investment  must  be  equal  to  the  budget.   This,  then, 
becomes  a  fixed  end  point  problem  as  described  in  reference 
(21). 

The  performance  equations  for  a  typical  interial  node 
as  shown  in  Fig.  5  can  be  written  as  follows: 

xn,m  m    (xn,m-l  +  xn-l,m  +  vn,m)bn,m  =  ^n.^.m         {6l) 
xn,m  m    (xn,m-l  +  ^-l.m  H  ^^m,  =  AIn.»(1..e».»j  (62) 


K"Am  L,n'mCt 


rn,m, n,m, 2 


vn,rr.   __n,mTn,m-   n,m  n,m   A12   1    t,   n,m  n,m. 


_i_   Kn  m 
Tn,m    ±3 
Ll 


+  x^"1  (63) 


v5 


Xg'm  -  K^nL^'mCtAIn'm(l-e^m) 

vn,rr;  Tn,m 


;2  ^2 

r. , ::. 

^2 .  „n,m 

Tn,r.  T  "23 
^2 


t  f.-rnm/i  -n'm\l  2   ,  vn,m-l       ,,.  , 

-  LAII''"(l-«3   )J    +  X6'  (o4) 
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-:...:  .  dn,m  4  er.,r.  +  xn,m-l  (65) 

where,   AIn'ffl  -  X^'01"1  +  X^1'1"  1    Vn>m  (66) 


and  vertical  links  respectively  at  node  (n,m). 

Since  tctal  investment  is  a  fixed  amount,  the  objective 
function  becomes: 

S  -  xfU   +  xJ'M  (67) 

The  Hamiltonian  function  is 

T,n, m    ,n,m  ..n,m  ,  „n,m  „n,m  ,  ,,n.m  vn,m   r.n,m  ,.n,m 

H       -   Z,      k-^  +  Z2      X,      +   Zj      lr      t  Z^      A£ 

+  Z^'r'  X2'm  (66) 

The  boundary  conditions  are  given  as  follows: 


.0,0  _  „o,o  _  ,-o,o  ^  .,0,0  _  vo,o 
cl   *  H       ~  A5       "  H       ~  "7       ~ 


,N,M  =  7N,M 


(69) 


(70) 


Z5'"  -  z^>r-  =  1  (71) 


XK,K-1  .  XN,M  =  GI  (72) 

Since  X~'M  is  fixed,  Zy'  '  becomes  an  unknown.   However, 
'with  xN,I''_1  fixed,  Zy'*"_1  can  be  obtained  by  solving  the 
following  equation: 


^f  - .   ICi Ct ^K,::-1)2  -  4,K_1  ■  ° (73) 

(.N,K-1  +  Kii"     ' 
Ll 

^N'M-1  =  GI  -  4-K-2 

Therefore,  zi^1  -  ^t_ _  (74) 

(7* +  kN.M"1) 

1   LN  M-l    A13 
Ll 

Since  Z"'"-1  =  ^n'm   =  Z*>m 

zn,rr.  _  ZM,M-1         for  ali  (n>m)  (?5) 


The  values  of  the  adjoint  variables  are  as  follows: 
Zx'    =  z2      zl  °3   T  Z2   ll-e3   ]  T  Kll  °3  Ll 


K^'m  i?,m  ct 

yn,mn    „n,mwn,mr  ,  -==■ = n,m,   n,ir.,2 

K2i    (1"93      )L2     Ct  T  2  e£>ffl  AI        (fa3      » 

Tn^E  ">    Klj 

Ll 
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..  Ln,ra  c 

4  2  _22 2 tA1nfB(1.en,B)2 

,r.  ,rr.  3  (7o) 


-2   n  Kn,m 


„n,m   „n,ir. 


j,a  -  Zg'™  =  1        for  all  (n,m)  (77) 


KK,M-lr  ,,TN,M-1, 


2 


Z">l:l  »   i2  >;  ,,^ for  all  (n,m)  (78) 

(  GI"X7'  "  +  rN.M-IjZ 

1    TN,K-1   rv13   ' 

Ll 


The  necessary  conditions  for  S  to  ce  a  local  minimum  is  that: 


3kLl!!  .  0  0  <en,m<  GI  _  xa,a-l  _  rf.m-1  (79) 

aan,m  5  6 


Oon  '  ^  °  -*- 


i^^         o<e^m<i  (8i) 


,     f  ,.n,rr.   n.rr.   .n.ffi  1  .      ,  . 
when  16,'  ,  8,'  ,  6,'   j  is  an  interior  point,  or 

h  '  =  minimum  (82) 

(    n  n"i   ^  in   ~  n ,  m^ 
when  Id''  ,  S_   ,  9,   is  at  a  boundary  point  of  the  con- 
straints. 


A  6 


equations  (61)  to  (66)  and  equation  (77) 
into  equation  (68)  and  taking  derivatives  with  respect  to 

various  decision  variables,  the  following  equations  are 
obtained : 


,.  n  ,  r.  / ,  T  r. ,  m  ,  n  .  r.  \  2r 
itin,m    K, i  (AI  '  o,'  I    C 


A  ,r.  ,rr.      n,m  7 

0ol       ,el      „n,m.  2 

(7^  +  K13  > 

Ll 


zV-1  (83) 


a0n,rn   "    ,n,m         -> 

lTn,ir,   4   K23  ' 
L2 


Z1^-1  (34) 


-  (Z?'m  -  z5'n)Aln»"  +  (K^X"m-K^'m)  CtAIn'm 


Kn,r,(A:r.,S)  28n,rr.Ln,m^ 
.  r. ,  bi 
-n,^    +   K13 


Kn,m(AIn,mj  2  (1.en,m)  Ln,m^ 
e2'r'  n  t 

*-  J.  T/Hj.L 

Tn,m         23 


(85) 


Taking  derivative  of  equation  (85)  with  respect  to  9"'m,  we 

obtain 
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,2Hn,m      K?-m(AIn'ra)2L?'niC  Kr  :n(AIn,m)2Ln,mc 

•2-i: 7  -  2  _i~ ± i  +  2  -22 ^ 1        (36) 

[dei'   )     61    .   „n,m  92     ,   Kn,m 

Tn,n      13  Tn,m      23 

1  L2 

Settinr  equations  (83)  and  (84)  equal  to  zero,  we  obtain  the 
following: 


n.m  ~*-*i?'   ct 

6   '    =    t     „r     -i    AI 


n.m  _n,m  Tn,ui   „n,m  Tn,m  /n,i 

o  '   L  '   -  K  '   L  '  (87) 


T^T"  Ki    °3   *1   -  "13  "l 


n,m 
n.m  _  K22  ct  ,Tn,m/n  ,n,m,Tn,m  vn>m   Tn,m  ,ai}1 

'?    ~   M  f  i   AI    (!-"-,   )L?    -  K9,   L9  (88) 


The  optimum  seeking  procedure  developed  for  this  problem  is 
as  follows : 

1.  Assume  a  set  of  decision  variables  (  8^>m,  e"'17"',  t^*111}  . 

2.  Calculate  values  of  Xn'm,  i  =  1,  2,  5,  6,  7  and  AIn'm 

i 

starting  at  n=m=l  and  proceeding  to  n=N,  m=M. 

3.  a.)   For  the  first  iteration,  calculate  Zy'   *~  by 

equation  (74)  with  the  above  X. '  '  and  AI  'm   values 
and  go  to  step  4. 
b.)   For  the  second  and  the  following  iterations,  cal- 
culate 2^'    by  equation  (74)  with  the  above  X?,m 
and  AIn,ir'  values.   This  z''    ~      value  is  then  com- 
pared with  the  value  obtained  in  the  previous 


iteration.   If  the  two  values  are  sufficiently 
close,  proceed  to  step  6.   If  they  are  not  suf- 
ficiently close,  proceed  to  step  4. 
Calculate  new  values  of  h^'^   and  82'   using  equations 
(87)  and  (38)  and  check  the  boundary  conditions. 
Return  to  step  2. 

With  the  above  9^'   and  X?,m  values,  calculate  Z^1'" 
and  Zp'^1  starting  at  n=N,  m=M  and  proceeding  backward 
to  n=m-l  by  the  use  of  equations  (70)  through  (76) 
Using  the  above  values  of  Xn'm  and  zn,m,  calculate 

— 4r-z:  arc  _  _  „  through  the  use  of  equations  (85) 

9\'  Ob,  '  )* 

and  (86). 

Adjust  the  values  of  8n,tn  by  adding  an  air.ount  equal  to 
A  ,  where 


ann 

,m 

3  6? 

,m 

a2:,:11 

,m 

A  =  - 

3' 

,  -  .  n  ,  m ,  2 
(3«3   ) 

and  check  the  boundary  conditions. 
9.    Return  to  step  2  and  repeat  the  procedure  until  the 
value  of  the  objective  function  [equation  (67))  is 
sufficiently  close  to  the  previous  value  to  indicate 
adequate  convergence. 
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In  the  case  where  a  minimum  level  of  service  is  to  be 
provided,  the  minimum  investment  can  be  treated  as  the  exist- 
ing facilities.   The  problem  can  then  be  solved  by  the  general 
method  without  chancing  the  algorithm.   In  other  words,  when 
the  values  of  K^j '  are  less  than  the  minimum  required  in- 
vestment ,  set  them  equal  to  the  minimum  investment  and  de- 
duct the  difference  from  the  total  budget. 

The  above  formulation  provides  solutions  to  a  single- 
quadrant  network,  single-copy  problem.   To  solve  a  multi- 
quadrant  network,  multi-copy  problem,  the  procedures  de- 
veloped by  Snell,  et .  al.  (23,  24)  can  be  employed.   How- 
ever, due  to  the  capacity  restriction  of  the  available 
computer,  this  extension  has  not  been  accomplished. 
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EXAMPLES  AND  DISCUSSION 

Several  examples  under  different  investment,  conditions 
are  presented  in  this  section  to  demonstrate  the  use  of  the 
model. 

A  hypothetical  network  was  developed  as  shown  in  Fig.  6. 
Node  U,i)  was  assumed  to  be  the  centroia  of  the  CBD.   Peak 
hour  trins  which  are  produced  in  the  other  zones  and  destined 
to  the  CBD  are  also  shown  in  the  figure.   All  links  have  an 
equal  length  of  one  mile.   The  area  was  divided  into  two 
parts  by  a  diagonal  line  which  passes  through  nodes  (1,A.) 
and  (/.,1).   The  lower  r>art  which  is  adjacent  to  the  CED  was 
assumed  lo  be  densely  developed.   The  upper  part  was  assumed 
to  be  less  densely  developed.   Assuming  the  maximum  speed  in 
the  densely  developed  area  to  be  60  mph  and  in  the  less 
densely  developed  area  70  mph,  minimum  travel  times  in  these 
two  areas  become  O.JI67  hour  per  mile  and  0.0143  hour  per 
mile  respectively.   Single  line  links  represent  existing 
local  streets  and  double  line  links  represent  existing 
arterial  streets. 

Innut  data  for  the  models  are  summarized  in  Table  1. 
Values  of  K.„  and  K   are  also  indicated  in  Fig.  7  and  Fig.  S 
respectively.   Since  construction  cost  and  right-of-way  cost 
will  not  be  the  same  in  each  area,  two  values  of  K.,  were 
assigned  to  the  links  even  though  these  links  represent  the 
same  type  of  facilities.   For  the  same  reason,  in  the  link 
investment  constraint  model,  links  have  different  values  for 
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2    3 
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incut  traffic  in  vph 

all  links  have  an  equal  length  of  1  mile 


(1,4) 


(2,4) 


(3,4) 


.(4,4) 
13,000 


Fig.  6  Hypothetical  Network  and  Peak  Hour  Traffic 
Distribution 
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Table  1   Input  Data  of  txample  Problems 


53 


. .  1  '0003 

0.00006 

0.00006 

-4 

u-\ 

^o 

3 

O 

o 

o 

H 

3 

• 

3 

3 

O 

3 

o 

o 

o 

o 

: 

o 

o 

0.00005 

o 

0.00006 

d 

0.00010 

3 

•JA 

u^ 

to 

W\ 

o 

3 

3 

t-\ 

o 

O 

O 

o 

- 

o 

O 

3 

o 

3 

O 

o 

3 

0 . 00006 

° 

o . oooos 

d 

0.00015 

3 

o 

O 

l/\ 

u-% 

o 

rH 

rH 

cv 

c 

3 

c 

3 

7 

3 

3 

3 

: 

3 

o 

3 

3 

o.oot 

3 

0. 00015 

3 

0.00020 

3 

'fig.    7  K.  „  Values  For  The  Example  Problems 


bk 


3 

r-i 

1C 

«0 

10 

15 

o 

r-i 

H 

12 

ITS 
H 

12 

X  5 

r-i 

Fig.  8  K.-,  Values  For  The  Example  Problems 


55 


maximum  and  minimum  investment  levels.   The  derivation  of 
these  data  is  discussed  in  Appendix  A.   Time  cost  (Ct)  is 
assumed  to  be  $1.55  per  hour  per  vehicle  as  suggested  by 
AASHO  (8). 

Example  1 :   Theoretical  optimal  system. 

Suppose  we  are  planning  for  a  completely  undeveloped 
area  where  no  facilities  exist  and  there  is  no  budget  limi- 
tation on  link  investment.   A  theoretical  optimal  system 
can  then  be  developed  to  accommodate  the  predicted  trip 
demand.   Using  the  formulation  of  "Investment  With  No 
Budget  Constraint"  and  letting  K"'m  =  0,  for  all  (n,m), 
the  resulting  system  is  shown  in  Fig.  9.   Notice  that  the 
system  forms  a  shortest  path  tree  in  which  only  one  route 
is  built  for  each  origin-destination  pair  and  all  trips  are 
assigned  to  this  route.   This  result  coincides  with  the 
analysis  discussed  in  page  33  which  shows  the  linear 
characteristic  of  the  problem  under  no  limit  condition. 

Example  2:   Optimal  investment  with  upper  and  lower  limits 
on  link  investment. 

The  hypothetical  network  shown  in  Fig.  6  is  to  be  im- 
proved with  the  following  conditions: 

1.  No  system  budget  limit. 

2.  A  minimum  level  of  service  (arterial  street)  is  to 
be  provided  for  the  entire  area. 
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3.    Roadway  space  obtainable  in  restricted. 

nvestmont  limits,  tj.min.  and  B.max.  ,  associated  with 
conditions  2  and  3  are  listed  in  Table  1.   The  formulation  of 
this  problem  has  been  aeveloped  in  the  previous  section  under 
the  category,  ''Investment  with  no  budget  constraint:   special 
case  III." 

The  results  are  shown  in  Fig.  10.   Note  that  with  the 
minimum  level  of  service  provided  for  the  entire  area,  trips 
are  assigned  rather  uniformly  to  take  advantage  of  all 
facilities.   When  traffic  is  focused  on  the  CBD,  a  space 
limitation  is  in  effect  which  forces  traffic  to  split  and 
enter  the  CBD  area  from  two  directions.   Considering  existing 
facilities  as  part  of  the  cost,  total  cost  becomes  §2,875.99 
(2,603.99  +  272.00).   Comparing  this  cost  with  the  total  cost 
in  example  1  ($2,319.86),  the  difference  is  only  about  two 
percent.   This  indicates  that  providing  a  minimum  level  of 
service  might  be  desirable  in  an  urban  area. 

Example  3 :   Investment  with  fixed  node  budget. 

This  is  the  fixed  node  investment  problem  as  formulated 
in  a  previous  section.   Investment  for  each  node  (SI  '  ')  is 
listed  in  Table  1.   Consider  the  area  as  completely  undeveloped 
K^m  =  0),  the  resulting  system  is  shown  in  Fig.  11.   Regard- 
less of  the  investment  level,  as  compared  with  previous  ex- 
amples, travel  Lime  costs  in  this  problem  are  greater  than 
in  the  previous  examples.   This  result  demonstrates  that 
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improper  allocation  of  funds  could  be  very  costly.   It  also 
indicates  the  advantage  of  area-wise  transportation  system 
development  which  is  carried  out  by  a  single  authority. 

.■:-;r.:-  ole  I.  :   Investment  with  fixed  system  budget. 

The  hypothetical  network  as  shown  in  Fig.  6  is  to  be 
improved  with  a  total  system  budget  of  $300  (GI  =  300, 
equivalent  peak  hour  budget).   The  resulting  traffic  assign- 
ment and  link  investments  are  shown  in  Fig.  12.   Most  in- 
vestment ar>oears  to  be  made  along  the  shortest  path  trees 
as  obtained  from  example  1  which  are  also  the  routes  that 
will  cost  "Che  least  to  improve.   This  is  logical  since  the 
budget  is  substantially  less  than  that  required  by  a 
theoretical  optimal  system  (compared  with  example  2). 
Comparing  the  costs  with  those  obtained  in  example  2,  it  is 
evident  that  although  investment  cost  decreases  more  than 
}C   percent,  total  cost  increases  only  1.4  percent.   This 
again  coints  out  the  advantage  of  area-wise  transportation 
system  development. 

In  order  to  verify  the  optimality  of  the  results,  two 
procedures  were  used. 

1.   Assuming  a  new  set  of  decision  variables  (»i'  ', 

^  m   n  m  \ 
bj   ,  8,'  )  to  start  with,  each  problem  was  solved  once 

more.   The  results  were  then  compared  with  the  previous  ones. 

No  significant  differences  between  the  two  solutions  of  each 

problem  were  observed.   Total  costs  are  summarized  in  Table  2. 
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stile  2 


Number  of  Iterations,  Approximate  Computing  Time  Used  And 
Total  Costs  for  Example  Problems 


Example  No. 

Starting  Point 

Number  of 
Iterations 

Time  Used 
(min. ) 

Total  Cost 
(?) 

1 

,,n.m    ,  „ 
B3   "  °'7 

15 

20 

2,819.36 

n ,  in   ,  „ 
ti-,   =  0.3 
3 

25 

30 

2 , 819 • 87 

2 

n ,  m 
9,   =0.7 

18 

25 

2,603.99 

,.n,m 

83   =  J .  3 

18 

25 

2,604.39 

3 

8  '  =  0.3 
i 

11 

15 

3,112.92 

n,m 
8-3   ■=  0 . 7 

12 

16 

3,112.95 

L 

cji   — -L?  bo   — > 
x-,r.  ,m   -~, 
93'  =  0. j 

16 

120 

2,639.38 

en,m=io  eg'^io 
^  e^-  0.7 

20 

160 

^,o_>7  -j° 

63 


Pinal  results,  in  the  original  computer  output  format,  are 
presented  in  Appendix  C.   Although  this  is  not  a  rigid  proof 
of  the  optimality,  it  does  indicate  that  the  results  are 
properly  converged  and,  therefore,  are  likely  the  optimal 
solutions . 

2.   Ore  decision  variable  was  selected  arbitrarily  and 
its  value  was  changed  from  one  percent  to  ten  percent  (some- 
what arbitrarily  but  related  to  its  original  value).   Keeping 
the  values  of  other  decision  variables  unchanged,  the  total 
cost  was  calculated.   This  cost  was  then  compared  with  the 
one  previously  obtained.   All  perturbations  resulted  in  a 
higher  cost  which  indicates  that  the  results  obtained  from 
this  model  are  at  least  very  close  to  a  local  minimum.   The 
costs  resulting  from  several  perturbations  for  each  problem 
are  shown  in  Appendix  C . 

Due  to  storage  limitations  of  the  available  computer 
(IBM  1620),  no  attempt  was  made  to  apply  this  technique  to 
a  more  complex  network.   The  above  examples  are  restricted 
to  one  quadrant,  single  copy  problems.   Therefore  only  a 
limited  comparison  of  the  results  to  the  real  world  is  pos- 
sible at  this  stage  of  development. 

Computer  programs,  one  for  each  budget  condition,  written 
for  use  on  the  IBM  1620  computer  are  presented  in  Appendix  3. 
The  number  of  iterations  and  approximate  computing  time  used 
for  each  example  problem  are  summarized  in  Table  2. 
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SIONS 

A  new  technique  for  the  analysis  of  transportation 
systeir.  investment  problems  has  been  presented  in  this  study. 
Considering  each  node  of   a  rectangular  network  as  a  stage, 
the  discrete  maximum  principle  was  utilized  to  formulate  a 
transportation  system  model.   Investment  models  under  dif- 
ferent investment  conditions  were  investigated  and  search 
procedures  were  developed  to  obtain  the  optimal  investment 
policy  and  to  assign  trips  to  the  network.   This  provides 
a  broad  application  of  the  model  to  solve  various  problems 
which  have  specific  investment  restrictions. 

As  opposed  to  linear  programming  models,  this  model  is 
capable  of   solving  transportation  system  investment  problems 
with  travel  time  being  non-linearly  related  to  traffic  volume 
and  investment  cost.   The  formulations  presented  in  this 
study  are  applicable  to  one  quadrant  network,  single  copy 
problems.   With  minor  modifications,  the  technique  would  be 
equally  applicable  to  a  more  complex  network  once  a  larger 
computer  is  available. 

The  optimum  seeking  procedures  appear'  quite  efficient 
and  yield  reasonable  results  as  shown  by  the  example  problems. 
Although  direct  comparison  of  the  results  with  the  real  world 
is  net  feasible  with  the  limitation  of  computer  capacity,  the 
model  does  represent  a  significant  step  toward  more  realistic 
analysis  of  transportation  systems. 
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Although  the  travel  time  function  derived  in  this  study 
could  be  further  improved  and  the  objective  function  might 
include  additional  variables,  this  study  has  demonstrated 
the  usefulness  and  the  ability  of  the  discrete  maximum 
principle  in  solving  this  type  of  non-linear  optimization 
problems.   It  also  indicates  that  the  discrete  maximum 
principle  could  be  a  powerful  tool  in  transportation  system 
analysis . 
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HECOKMISNDATIOKS  FOR  FURTHER  RESEARCH 

The  relationships  between  travel  time,  trafi'ic  volume 
and  investment  cost  are  complex.   Although  the  non-linear 
travel  time  equation  developed  in  this  study  is  considered 
to  be  more  realistic  than  a  linear  approximation,  there  is 
undoubtedly  room  for  further  improvement.   Two  improvements 
which  might  be  considered  are:   (a)  the  desirability  of 
capacity  restriction  on  links,  (b)  the  relationships  of 
free  flow  travel  time  to  investment. 

Previously,  it  was  mentioned  that  an  approximate  pro- 
cedure has  been  developed  to  solve  multicopy  problems. 
However,  this  procedure  will  become  invalid  when  the  system 
budget  is  restricted,  which  is  not  uncommon  in  the  real 
world.   Therefore,  to  develop  a  useful  model,  an  improved 
technique  is  required.   One  possible  approach  is  to  con- 
sider each  copy  as  a  large  stage  and  each  node  becomes  a 
small  stage  inside  the  larger  stage.   The  objective  function 
of  each  copy  becomes  the  Hamiltonian  of  the  larger  stage. 
By  this  concent,  the  discrete  maximum  principle  could  be 
utilized  to  develop  a  complete  model  for  multicopy  tran- 
sportation problems. 
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APPENDIX  A 

VALUES  OF  CONSTANTS  IN  UNIT  TRAVEL  TIKE  EQUATION 

Unit  travel  time  has  been  expressed  as: 

..    K2 
t  ■=  K,  +  


t  =  unit  travel  time  (hr/mi/veh) 
Kj_  =  free  flow  travel  time  (hr/mi/veh).   The  magnitude 

depends  on  the  maximum  speed  obtainable  or  regulated. 
K-  =  coefficient  of  improvement  ( dollar-hr/mi'yveh  )  • 

Its  magnitude  depends  on  link  location  and  reflects 

the  difficulty  of  improvement. 
K,  •=  existing  investment  (dollar/mi/hr ) 
o  =  equivalent  hourly  investment  per  unit  length 

(dollar/mi/hr) 
V  -  traffic  volume  per  unit  time  (veh/hr) 

In  this  section,  a  set  of  K  values  are  derived  from  real 
world  data  reported  by  other  researchers.   The  purpose  of 
this  section  is  two  fold: 

1.  To  justify  the  fitness  of  the  equation. 

2.  To  obtain  a  set  of  K  values  for  the  example  problems. 

Values  of  K-, 


The  Kt  value  is  equal  to  the  reciprocal  of  the  maximum 
speed  obtainable  or  regulated  in  each  area.   Several  common 


e  shown  in  Tabic  A-l. 


Table  A-l 


mum  Speed           K]_  Values 
(mph)               (hr/ir.ile) 

70 

0.0143 

60 

0.0167 

50 

0 . 0200 

For  the  example  problems,  maximum  speeds  were  assumed 
-„o  te  70  r.ph  in  less  densely  developed  areas  and  60  mph  in 
densely  developed  areas.   The  K-j_  values  are  therefore, 
0.0113  hours  per  mile  and  0.0167  hours  per  mile  respectively. 

Values  of  K9 


1.   Near  CBD  Area: 

The  average  cost  of  an  8  lane  freeway  near  the  CBD,  as 
estimated  by  Aitken  (19),  is  $15,500,000.00  per  mile.   As- 
suming 30  year  life  and  6%   interest,  annual  cost  is  equal 
to  $1,130,000.00  per  mile.   If  we  further  assume  peak  hour 
traffic  is  10$  of  daily  traffic,  the  equivalent  peak  hour 
cost  becomes: 


$1,130,000  X 


1 


— ] —  =  v3l4  ser  mile  cer  hour 
10 


This  freeway  can  handle  1100  vph  per  lane  at  unit  travel  time 
of  0.02  hr/mile.   Assuming  Kx  =  0.0143  hr/mi/veh  (70  mph 


Spec  .        ws: 

:  ...:..;■  h   "■•    (1100  X  8)  -  0.020 
3H 


K,  -  0.00207  dollar-hr/mi2/veh2  (A-2) 

Using  Haikalis'  data  and  adjusting  for  the  downtown  area 
(19),  an  arterial  street  with  2,000  vph  volume  at  unit  travel 
time  of  O.G333  hour  per  mile  costs  $3,400,000  per  mile  or 
$250,000  per  mile  annually.   Equivalent  peak  hour  cost  be- 
comes : 

$250,000  X  — t —  X  — —  -  $69-5  per  mile  per  hour. 
3o0     lo 

Assuming  K]  ™  0.025  hr/mi/veh  (40  mph  speed),  K  is  derived 
as  follows: 

0.025  +  — — —  2,000  =  0.0333, 
69.5 

K,  =  0.0'8O2S£  dcllar-hr/mi2/veh2  (A-3  ) 

2.   Average  Urban  Area: 

The  overall  average  cost  for  an  8  lane  urban  freeway  is 
$5,001  ,000  per  mile  as  estimated  by  Koskowitz  (26).   Assuming 
30  year  life  ar.d  6%   interest,  equivalent  peak  hour  cost  be- 
comes : 

$5,000,000  X  0.0726  X  — = —  X  — L-  =  $101  ser  mile  oer  hour. 
360    10 
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.  ■..;■.,■'  in     "Hi  :hway  Capacity  Manual"  (25) ,  a 
typical  freeway  will:  '."  tup  average  highway  speed  cor,  handle 
1800  vph  per  lane  at  a  speed  of  1,5   mph .   The  K2  value  is 
derived  as  follows: 

K-,  =  0.0143  hr/n.i/veh 

0.0143  +  Ji_   (1600  X  8)  =  0.0222 
101 

K2  =  0.0000553   doliar-hr/mi2/veh2 

Summarized  from  "Automobile  Transportation  Systems: 
Cost  Characteristics"  (27),  Table  A-2  shows  relationships 
amor.?;  volume,  average  speed  and  cost  for  three  types  of 
urban  roads.   Using  these  values  and  the  assumed  maximum 
speeds  arid  average  lanes,  Table  A-3  is  obtained.   The  K2 
values  are,  then,  derived  as  follows: 

K2 

local  street:   0.0286  +  — —  X  1000  =  0.05 

K2  =  0.000227  dollar-hr/mi2/veh2       (A-5) 

K, 

arterial  street:   0.025  +  _s: X  28u0  =  0.308 

19.2 

K2  =  0.0000398  dollar-hr/mi2/veh2      (A-6) 

K, 
freeway:   0.0143  +  -4 —  10S00  -  0.0182 
161 


K2  -  0.0000582  doilar-hr/mi2/veh2      (A-7) 


76 


Table  A-2  Cost   Characteristics  of  Urban  Highways 


Practical  Capacity 
(vph/lane ) 

Average  Speed 
(mph) 

ROW  Cost 
($/raile) 

Construction  Cost 
(t'/miie  ) 

Total  Cost 
($/mile) 


Local  Street  Arterial  Freeway 

5U0       700  1800 

20      25-40  45-^5 

250, 00U  450,000  4-8  million 

300,000  500,000  4-6  million 

550,000  950,000  8-14  million 

source  =  Ref.  (28) 


fable  A-3  Cost  and  Travel  Time  of  Urban  Highways 


No.  of  lanes 

Total  Volume 
(vph) 

Average  Speed 
(mph) 

Total  Cost 
(S/mile) 

Equivalent  Peak  -hour 
cost   (S/mile) 

Assumed  maximum 
speed    (mph) 

Minimum  unit 

travel  time  (hr/mile) 

Average  Travel 
Time  (hr/mile) 


cal  Street 

Arterial 

Fr 

eeway 

2 

4 

6 

1,000 

2,800 

10,800 

20 

32.5 

55 

550,000 

950,000 

8-14 

million 

11.1 

19.2 

161-282 

35 

40 

70 

0.0286 

0.025 

0.0143 

0.05 

0.0308 

0.0182 

source  =  Ref .  (28) 
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0.0143  -i  -   luSoo  -  0.0182 
K2  «=  0.000102  dollar-hr/mi2/veh2 
Rural  Area: 


(A- 


Cost  aata  for  rural  highways  is  not  generally  available. 
However,  the  cost  of  a  rural  freeway  nay  be  assumed  as  equal 
to  the  lowest  cost  of  a  freeway  in  an  urban  area. 

Or.  this  basis  an  8  lane  freeway  will  cost  about  ;;;3,00u,000 
per  mile  (27).   Using  Fig.  3-38  in  the  "Highway  Capacity 
Manual"  (25),  a  typical  freeway  with  70  mph  average  highway 
speed  can  handle  1800  vph  per  lane  at  45  mph  speed.   Equiv- 
alent peak  hour  cost  becomes: 


$3,000,000  X  0.0726  X  ri-  X  ~ 
3o0    10 


$60.5  per  mile  per  hour 


The  K-  value  is  derived  as  follows: 

K2 
0.0143  +  7— —  (1800  X  8)  =  0.0222 
60.  5 


K, 


.00003322  dollar-hr/mi^/veh' 


(A-9) 


Excluding  equation  A-5,  Kz  values  are  summarized  in  Table 
A-4. 

Table  A-4 


type  01  Ares 


Range  of  K5  Value 


C3J 


0.000207-0.000288 


Average   Urban  Ar< 


Rural 


0 . 0000398-0 . 000102 

0.000^332 
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..,..'.    .... 
free-way  costs  as  shown  in  Table 
airly  consistent  in  each  area. 
3d  correlation  between  the  equation 
world  situation. 

The  ':. ,    value  represents  the  existing  facilities  in  terms 
of  cost  oar  mile  per  hour.   Equivalent  peak  hour  cost,  for 
each  type  of  road,  derived  in  the  previous  sections  indicates 
the  avera  ; ;  ral\  '  .  . 

The  K  val  es  usee  in  the  example  problems  are  su.-r.rr.arized 
Fable  1. 


APPENDIX   B 
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RESULTS  OF  PERTURBATION  ,  PROBLEM  1 

Decision  Variables 
being  Perturbed 

Original 
Vaiues 

Values  after 
Perturbation 

Resulting 
Total  Costs 

^ 

1 . OOOO 

•  9500 

622818.51 

el)2 

44.0100 

40.0000 

2820.25 

92,3 

.0000 

.1000 

2819-95 

#* 

15.2400 

16.0000 

2819-88 

.}.» 

1.0000 

.9500 

31390309.00 

fi3,4 

92 

19.6300 

18.0000 

2820.01 

e^-1 

11.1300 

12.0000 

2819.91 

e^3 

211.2000 

215.0000 

2819.91 

ORIGINAL  TOTAL 

COST  =  $2,819 

.86 

„ 

RESULTS  OF  PERTURBATION  ,  PROBLEM  2 

Decision  Variables 
Being  Perturbed 

.1   T 

Original 
Values 

Values  after 
Perturbation 

Resulting 
Total  Costs 

H1-1 

V 

2.0000 

2.2000 

2604.10 

bJ.* 

•  3834 

.4000 

2604.01 

e2,2 

24.0300 

23.0000 

2604.03 

e2,3 

•  3743 

.3600 

2604.05 

93.1 

2.0000 

2.2000 

2604.06 

e3,3 
3 

•  4307 

.4200 

2604.07 

e^.i 

6.0800 

6.5000 

2604.00 

a*. 3 

X 

85.0000 

84.0000 

2604.62 

ORIGINAL 

TOTAL  COST  =  $2,603 

.99 
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RESULTS  OF 

PERTURBATION  ,  PROBLEM  3 

Decision  Va 

riables   Original 

Values  after 

Resulting 

Being  Perturbed       Values 

Pe; 

rturbation 

Total  Costs 

«$.! 

.5231 

.5000 

3112.95 

e£.3 

21.7300 

22.5000 

3112.91 

e2,2 

2 

31.8000 

31.0000 

3112.92 

e2,3 
3 

•  3174 

.3100 

3112.95 

93,1 
1 

29.6600 

31.0000 

3112.92 

e3,3 
3 

.4083 

.4000 

3112.96 

ORIGINAL  TOTAL  COST  =  $3,112 

.90 

Decision  Variables 
Being  Perturbed 

Original 
Values 

Values  after 
Perturbation 

Resulting 
Total  Costs 

V 

.5721 

.5500 

V 

.5721 

.5500 

2639.49 

e^.2 

8.0040 

.8800 

92-1 

5.2560 

4. 4560 

2654.22 

°3'2 

.2672 

.2800 

,2,3 

.2672 

.2800 

2639.52 

eg.* 

4.7630 

5 . 0000 

e3,i 

15.2200 

14.9800 

2639.46 

93,1 

.4999 

.4600 

e^.i 

.4999 

.4600 

2639.65 

8*.  2 

.9788 

1.0000 

eA,2 

23-3100 

23.2900 

2639.44 

9U,3 

23.3100 

23.0000 

108 . 9000 

109.2100 

2639.45 

ORIGINAL  TOTAL 

COST  =  $2,639.38 
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LIST  OF  SYMBOLS 

Discrete  Maximum  Principle 

X  state  variable. 

9  decision  variable. 

T  transformation  operator. 

n  the  n-th  stage. 

N  the  N-th  stage  or  the  total  number  of  stages. 

s  total  number  of  state  variables  in  each  stage. 

r  total  number  of  decision  variables  in  each  stage. 

H  the  Hamiltonian. 

z  adjoint  variable. 

c  constant  in  objective  function. 

S  objective  function. 

General  Formulation  of  the  Problem 

Xn,m     state  variables  representing  flows  from  node  (n,m). 


-.rum 
3j 


decision  variables  representing  investments  on  links 
leaving  node  (n,m). 
K-;{      free  flow  travel  time  on  links  leaving  node  (n,m). 

K'      coefficient  of  investment  on  links  leaving  node  (n,m), 

K.'i  '     existing  investment  on  links  leaving  node  (n,m). 

L.'      link  length  on  links  leaving  node  (n,m). 
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t  .'  unit  travel  time  on  links  leaving  node  (n,m). 

where,  j  =  1,  for  horizontal  link. 

j  =  2,  for  vertical  link. 
n,m 
X3       state  variable  representing  the  total  investment  on 

horizontal  links  from  node  (1,1)  through  node  (n,m), 

X.  '      state  variable  representing  the  total  investment  on 


,n,m 


vertical  links  from  node  (1,1)  through  (n,m). 
state  variable  representing  the  total  travel  time 
cost  on  horizontal  links  from  node  (1,1)  through 
node  (n,m). 
X^'      state  variable  representing  the  total  travel  time 

cost  on  vertical  links  from  node  (1,1)  through  (n,m), 

„  n ,  m  .  ,  -,  .  ., 

X7       state  vanaDle  representing  the  total  investment  on 

both  links  from  node  (1,1)  through  node  (n,m). 
60'      decision  variable  representing  the  fraction  of  the 

vehicles  departing  node  (n,m)  on  the  horizontal 

link. 
C^       time  cost. 

Hn'm     Hamiltoniam  function  at  node  (n,m). 
V  'm     input  trips  at  node  (n,m). 
GI       total  system  budget. 
SIn,m    section  budget  at  node  (n,m). 
Z,'  ,  Z  '    ,  Z  '   adjoint  variables  associated  with 


vn,m  vn,m      vn,m       . .   , 
1   »   ?   *  * " ' '  »  7   respectively. 


7 

objective  function. 
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Al  STRACT 

Following  an  introduction  to  the  purpose  of  the  economic 
analysis  of  transportation  networks,  the  objective  and  cri- 
terion problems  arc  briefly  discussed.   A  single  objective  is 
selected  for  this  study.   A  review  of  the  literature  shows 
the  historical  development  of  the  methods  of  economic  analysis. 
Major  drawbacks  of  the  traditional  methods  developed  in  the 
past  are  described.   The  merits  and  demerits  of  mathematical 
models  are  presented. 

The  discrete  maximum  principle  is  introduced  with  a 
brief  review  of  its  recent  applications  to  transportation 
systems  analysis.   Starting  with  a  description  of  the  rela- 
tionship between  travel  time,  traffic  volume  and  invest- 
ment cost,  a  non-linear  total  travel  time  equation  is 
developed  which  expresses  travel  time  as  a  function  of 
traffic  volume  and  investment. 

The  purpose  of  this  research  was  to  formulate  an  optimal 
network  improvement  model  (in  equation  forms  as  described)  by 
the  discrete  maximum  principle.   Utilizing  these  equations, 
optimum  seeking  procedures  were  then  developed.   Three  in- 
vestment conditions  were  considered  which  resulted  in  three 
different  sets  of  equations  and  two  slightly  different  ways 
of  seeking  the  optimum.   Three  special  cases  which  implied 
limits  on  link  investment  were  also  described.   This  formu- 
lation provided  a  broad  application  of  the  technique  to 


problems  with  various  constraints  and  assumed  conditions. 
Finally,  four  examples  were  presented  to  demonstrate 
the  usefulness  of  the  technique  in  different  investment 
conditions.   Derivation  of  the  data  used  in  example  problems 
and  the  computer  programs  developed  to  solve  the  problem 
were  presented  in  Appendix  B. 


